Stability in the evolution of random networks 



L.P. ChiQ C.B. Yang, and X. Cai 
Institute of Particle Physics, Hua-Zhong (Central China) Normal University, Wuhan 430079, P.R. China 

(Dated: February 2, 2008) 

With a simple model, we study the evolution of random networks under attack and reconstruction. 
We introduce a new quality, invulnerability I(s), to describe the stability of the system. We find that 
the network can evolve to a stationary state. The stationary value 7 C has a power-law dependence 
on the initial average degree (k), with the slope is about —1.485. In the stationary state, the degree 
distribution is a normal distribution, rather than a typical Poisson distribution for general random 
graphs. The clustering coefficient in the stationary state is much larger than that in the initial state. 
The stability of the network depends only on the initial average degree (k) , which increases rapidly 
with the decrease of (k). 

PACS numbers: 89.75.Hc, 87.23.Kg, 89.75.Fb 



I. INTRODUCTION 

Many systems can be represented by networks, a set 
of nodes joined together by links indicating interactions. 
Social networks [J, the Internet Q , food webs Q , trans- 
portation networks 0, 0, E| , and linguistic networks Q 
are just some examples of such systems. The investi- 
gation of complex networks was initiated by Erdos and 
Renyi H • They proposed and studied one of the earliest 
theoretical models of a network, the random graph. In 
a random graph, N labelled nodes are connected by n 
edges, which are chosen randomly from the N(N — l)/2 
possible edges. It is trivial to show that the connection 
probability is p = n/[N(N — l)/2]. The number k of 
edges connecting one node to others is called the de- 
gree of that node. The average degree of the graph is 
(k) =2n/N = p(N - 1) ~ pN if p < 1. The degree dis- 
tribution for a random network is given by a Poissonian 
distribution. 

Recently the increasing accessibility of databases of 
real networks and the availability of powerful comput- 
ers have made possible a series of empirical studies on 
complex networks. Thus, other two main streams of 
topics were proposed and investigated in depth. One is 
the small-world networks introduced by Watts and Stro- 
gatz [{J |l0( . Such networks are highly clustered like reg- 
ular lattices, yet have small characteristic path lengths 
like random graphs. The other is the scale- free networks 
proposed by Barabasi and Albert 0, 0, , based on 
two generic mechanisms, growth and preferential attach- 
ment. Those networks have scale-free power-law degree 
distributions. 

With increased threats of hacker attacks and routers 
malfunction, etc., research in the field of network ro- 
bustness has attracted much attention ^4|. Albert and 
his collaborators have shown that scale-free networks, at 
variance with random graphs and small-world networks, 
are almost unaffected by errors while vulnerable to at- 
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tacks. That is, the ability of their nodes to communi- 
cate is almost unaffected by the failure of some randomly 
chosen nodes, but the removal of a few most connected 
nodes can damage the networks. Recently the network 
efficiency of errors and attacks on scale- free networks has 
also been studied ^5|. In those studies, the damaged 
nodes are removed from the network. Thus the size of 
the network decreases with the evolution of the system. 

Based on the physics of network tolerance, using a sim- 
ple model on evolving network, we are trying to study the 
robustness of a dynamical evolving network. We will con- 
sider the reconstruction of the links of the damaged nodes 
instead removing the nodes from the network. Thus the 
size of the system remains unchanged. This is more closer 
to the evolution of real networks. The paper is organized 
as follows. The evolution model is presented in next sec- 
tion. Section III is devoted to analyze the results of the 
model, along with a definition of how we describe the 
stability of the network. In last section we conclude and 
give a brief discussion. 



II. AN EVOLUTION MODEL 

Despite the fine work of studies on network tolerance, 
little effort has been made on the reconstruction of the 
attacked network. The aim of our model is to investi- 
gate the stability of the network during the evolution in 
terms of attacks and reconstructions. In other words, the 
nodes damaged will not be removed from the network, in- 
stead they will be reconnected in certain way. Assume 
that all information on the former links of the damaged 
nodes has been lost, the damaged nodes have to be con- 
nected randomly to other nodes in the network again. In 
this way, we keep the size of the system constant. We 
do not try to consider a case with increasing size in the 
evolution, since the process of system's size increase can 
be much slower than the frequently happening damage 
and reconstruction. We try to investigate the effect of 
such a reconstruction on the evolution of the network. 
For this purpose, we first setup a random netweok with 
N nodes and connection probability p, then let the net- 
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FIG. 1: Plots of k max versus step s with N = 1000, p re 
p = 0.02. 
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FIG. 2: Plots of invulnerability I(s) versus s for network with 
N = 1000 nodes and connection probability p re = p — 0.02. 



work evolve according to the following rules: (i) Find the 
node with the highest degree k max . Since this node has 
highest connections to other nodes, it is most likely at- 
tacked and is the most vulnerable site in the network. 
For simplicity, we assume that only a node with high- 
est connections suffers attack and to be reconnected. If 
several nodes happen to have the same highest degree of 
connection, only one (randomly chosen) of them is as- 
sumed to be damaged in the attack, (ii) Reconnect this 
node with the other nodes in the network with the re- 
construction probability p re — p. Steps (i) and (ii) are 
repeated to a prefixed number of times. This number 
should be large enough to enable the system to reach 
(possibly existing) stationary state. In our simulation, 
it is chosen to be 10 million. The reconnection process 
represents the effort of reconstructing the network. We 
set the reconnection probability p re the same as the ini- 
tial construction probability p to reduce the number of 
parameters in the model. This set also implies that the 
information does not increase from initial construction to 
later reconnection. Thus, the evolution of the network is 
in fact a process of bing damaged and subsequent recon- 
structing, similar to the evolution of real networks. 

The interest in our model investigation is to find out 
whether there exists a stationary state in the evolution 
of the network, and if it exists, what are the properties 
of the stationary state and their dependences on the pa- 
rameters in the model. 



III. SIMULATION RESULTS 

The model defined in Sect. II has several nontrivial 
consequences. It can be easily seen that the network has 
a tendency to decrease the maximum number of connec- 



tions among the nodes at a long time scale, because the 
nodes with highest connections will be damaged and re- 
connected randomly. From the evolution rules, the total 
number of connections in the network will also decrease 
generally before a stationary state is reached. 

To get some idea on the safety of the network, we can 
have a look at the behavior of the maximum degree of 
connection of the network. In Fig. 1, we give a snap- 
shot of the maximum degree k max versus time step s for 
a network with N = 1000 nodes and connection proba- 
bility p re = p = 0.02. Apart from some fluctuations, it 
decreases in the evolution. From intuition, a node with 
less links to others will be attacked less frequently. Thus 
a network with smaller maximum connection degree is 
safer. To describe the safety of the network, we introduce 
a new quantity, the invulnerability I(s), which is analo- 
gous in definition to the gap G(s) in the Bak-Sneppen 
model fllilr7| . Considering an evolution of network with 
maximum degrees fc max (l), fc max (2), • ■ • , fc max (s), invul- 
nerability 7(s) at time s is defined 



J(s) = 1/Min{k max (i)} for i < s 



(1) 



i.e., the inverse of the minimum of all the maximum de- 
grees in the evolution before moment s. Initial value of 
I(s) is equal to l/k max (l). I(s) reflects the attack toler- 
ance of the network. When I(s) is small, the network is 
vulnerable to attack. Obviously from definition, I(s) is 
a non-decreasing function of evolution, and some fluctu- 
ations in fc max have been filtered out. 

Fig. 2 shows I(s) versus step s with network size 
N = 1000 and probability p re = p = 0.02. We observe 
that I(s) increases very quickly at small s but slowly at 
large s and finally reaches a constant value when s is large 
enough. The increase of I(s) indicates that the system 
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FIG. 3: The stationary value of I c as a function of network 
size N when (fc) = 20 fixed. 



FIG. 4: ln/ c versus ln(fc). A linear fit is shown as dash line 
with slope —1.485. 



is getting more and more safe in general by experienc- 
ing continuous attacks and reconstructions. This figure 
is similar to the envelope function of Bak-Sneppen evolu- 
tion model. Without interference from outside world, the 
network evolves to a stationary state. And the process 
takes place over a very long transient period. In Fig. 3, 
we present the stationary value I c as a function of the net- 
work size N under a fixed initial average degree (fc) = 20. 
We find I c stays unchanged at 0.2 with N shifted from 
100 to 5000. This result is interesting, because it shows 
that the stationary value I c depends not on the network 
size N and probability p re — p separately but through 
the initial average degree (k) = pN. Such a dependence 
can be expected considering the facts that the initial de- 
gree distribution has an average (k) = pN with variance 
a 2 = (k) = pN and that the average number of links to 
the reconstructed nodes is also (k) — pN. As a result, 
the properties of the evolution are mainly determined by 
the value of initial average degree (k) — pN. To get the 
relationship between I c and (k) we show in Fig. 4 I c as 
a function of (k) in a log- log plot. We find that I c has a 
power-law dependence on the average degree (fc), 
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where the exponent r is about 1.485. Fig. 4 illustrates 
that after the network has relaxed to the stationary state, 
the stability of the network will increase rapidly with 
the decrease of average degree (k) in the initial state. 
Thus, when the initial average degree (k) is small, i.e., 
less communications and interactions in the network, the 
system will be more stable. We would like to mention 
that the value of r = 1.485 results from the nonlinear 
interactions among the nodes. 

In order to offer a further information of the network 



in the stationary state, one can compare some structural 
properties of the network in the initial state with those 
in the stationary state. 

Because of the reconnection of a damaged node to 
other randomly chosen ones, the nodes with modest de- 
gree of connection may increase their links. As a result 
of the reconstruction, the distribution of connections k 
changes in the evolution of the system. To see how dra- 
matically the change is, we compare the distributions 
at initial state and the state in which I(s) remains un- 
changed. For this purpose, we construct a random graph 
with A^ = 10000 nodes and p = 0.005 and plot its degree 
distribution in the initial state by circles in Fig. 5(a). 
Because of the randomness of the initial connections, the 
initial degree distribution of a random graph is a Poisson 
distribution for large N [13J, 



P(k) 



-<fe) 



(k) k /k\ 



(3) 



with peak at about (A;) = 50. This theoretical expecta- 
tion is shown in the figure as the solid line. Then let the 
network evolve with the given simulation rules. When the 
system reaches the stationary state, we draw its degree 
distribution, as represented by squares in Fig. 5(a). We 
fit the distribution to a normal distribution, and find that 
the average degree is (k) = 17 and standard deviation is 
cr = 7. One can see that the mean degree of connection 
decreases by a factor 3 in the evolution. One may notice 
that a Poissonian distribution can be approximated by 
a normal distribution when the mean value (k) is large. 
For Poissonian distribution, the width of the distribution 
is determined also by the mean value {k) as a = \J (k). In 
a normal distribution, the width has no correlation with 
its mean value. Therefore, the stationary degree distribu- 
tion can not be described by a Poissonian distribution. 
To investigate the dependence of such a change on the 
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initial (fc), we do the same investigation for a network 
with N = 10000 nodes and p = 0.015. This time, the 
initial (k) is 150. The initial and final state distributions 
of connection are shown in Fig. 5(b), represented by cir- 
cles and squares, respectively. Both distributions can be 
fitted very well by Gaussian distributions, with (k) = 150 
for the initial distribution as expected and 77 for the final 
state. The width of the final state connection distribu- 
tion is now (7 — 18.5. From these two evolutions one can 
conclude that the change of connection distributions de- 
pends on the initial mean degree nontrivially. The higher 
the initial mean degree, the wider of the distribution for 
the stationary state. In fact, the two average degrees 
above in the stationary states satisfy the relation in Eq. 
(J2J with the same r for I c . One more observation is that 
the relative width a / (k) of the distributions increases by 
a factor of about 3 in the evolution for different initial 
average degree (k) . Because of the constant factor in the 
change of relative width in the evolution, the width for 
the stationary state can be smaller or larger than that in 
the initial state, depending on the value of initial average 
degree. 

Another feature with the evolution of the network is 
the emergence of isolated nodes. When the only one link 
a node has is a connection to the node with highest de- 
gree under attack, the node may have no connection to 
the network and becomes isolated after reconstruction. 
Needless to say, the number of isolated nodes depends 
on the evolution stage and the value of initial average 
degree (k). The larger the initial average degree (k), the 
less the number of isolated nodes at fixed steps since the 
evolution. Our simulation shows that the probabilities 
for a node to be isolated in the stationary state are 0.03% 
and 25% with initial degree (k) = 50 and 15, correspond- 
ingly. When initial mean degree (k) is small, many nodes 
are isolated, and the isolation of many nodes makes the 
nodes more independent and the system more stable. 

An important property of a network is its clustering 
coefficient. To know the behavior of the clustering coeffi- 
cient in the evolution, we need to calculate the clustering 
coefficient of the network in the initial state and in the 
stationary state. Let us focus first on a selected node i of 
the network. Suppose the node have fcj edges connecting 
its so-called nearest neighbors. The maximum possible 
edges among fcj nearest neighbors is ki(ki — l)/2. We 
use rii to denote the number of edges that actually ex- 
ist among those neighbors. The clustering coefficient of 
node i is defined as: 

° l = hih-l)/* (4) 

The clustering coefficient of the entire network is de- 
fined as the average over the whole network, 

C = ^Y. C >- (5) 

i 

We find that the clustering coefficient in the station- 
ary state is much larger than that in the initial state. 




FIG. 5: A comparison of degree distributions of the network 
in the initial state (circles) with that in the stationary state 
(squares) for (a) N = 10000 and p = 0.005; (b) N = 10000 
and p = 0.015. Solid lines are Gaussian fits. 

For a network with N = 10000 nodes, the clustering co- 
efficients in the initial state are equal to 0.005 and 0.015 
when the initial mean degrees (k) are 50 and 150, respec- 
tively, while the clustering coefficients in the stationary 
state are 0.23 and 0.1, respectively. 

From the fact that the network in the stationary state 
has the large clustering coefficients together with some 
isolated nodes and low average degrees, one can conclude 
that the system is driven in its evolution to a state com- 
posed of quite a few highly clustered small clusters. This 
is in sharp contrast with the network in the initial state 
when the system has a small clustering coefficient, almost 
no isolated nodes, and a high average degree. 

IV. DISCUSSIONS AND CONCLUSIONS 

In this paper, we study the evolution of random net- 
works under continuous attacks and subsequent recon- 
structions. We introduce a new quality, invulnerability 
I(s), to describe the safety of the network. A stationary 
state with fixed I(s) is observed during the evolution of 
the network. The stationary value of invulnerability I c 
is found to be independent of the network size N and 
the probability p re = p when the initial average degree 
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is fixed. I c shows a power-law dependence on the initial 
mean degree (k), with the exponent is about —1.485. 

We give further information on the evolution of the 
properties of the network. The first property is the evolu- 
tion of the degree distribution. The degree distributions 
of a network in both the initial and the stationary states 
are found to be normal distributions. After evolution, 
the peak position in the distribution shifts to lower con- 
nection degree while the relative width o~/{k) increases 
considerably. In the stationary state, the edges and de- 
grees in the whole network decrease a lot and quite a 
few isolated nodes appear. The second is the clustering 
coefficient. The clustering coefficient of a network in the 
stationary state is much larger than that in the initial 
state. 

In summary, the stability of the network is found to 
be related closely to the initial average degree (k) and 
has little correlation with the network size N and prob- 
ability p re = p separately. The stability of the network 
will increase rapidly with the decrease of initial average 
degree (k). The reason is that when (k) is small, the 
edges and degrees decrease a lot and more isolated nodes 
appear, which weaken the ability of nodes to communi- 
cate with each other and make the network more stable. 



From the isolated nodes and high clustering coefficient, 
we conclude that the network in the stationary state is 
composed of some highly clustered small clusters. 

Still, there are many issues to be addressed, such as the 
correlations and the fluctuations in k ma x in the evolution, 
especially after the stationary state. These fluctuations 
may tell us more the nature of the stationary state. The 
behavior of the average degree in the evolution also may 
shed some light on the dynamics of network. In addition, 
it is worthwhile to investigate whether the evolution of 
random networks demonstrate self-organized criticality 
(SOC), according to the similarity between the evolution 
of invulnerability J(s) with that of the envelope function 
G{s) in Bak-Sneppen model. All these topics can not be 
covered in this paper, and will be discussed later. 



Acknowledgments 

This work was supported in part by the National 
Natural Science Foundation of China under grant No. 
70271067 and by the Ministry of Education of China un- 
der grant No. 03113. 



[1] S. Wasserman and K. Faust, Social Network Analy- 
sis, (Cambridge University Press, Cambridge, England, 
1994). 

[2] R. Albert, H. Jeong, and A.-L. Barabasi, Nature 401, 
130 (1999). 

[3] R. J. Williams and N. D. Martinez, Nature 404, 180 
(2000). 

[4] W. Li and X. Cai, Phys. Rev. E 69, 046106 (2004). 
[5] L. P. Chi, R. Wang, H. Su, et al, Chin. Phys. Lett. 20, 
1393 (2003). 

[6] V. Latora and M. Marchiori, Physica A 314, 109 (2002). 
[7] R. F. Cancho and R. V. Sole, Proc. Natl. Acad. Sci. 100, 
788 (2003). 

[8] P. Erdos and A. Renyi, Publ. Math. Inst. Hung. Acad. 

Sci. 5, 17 (1960). 
[9] D. J. Watts and S. H. Strogatz, Nature 393, 440 (1998). 



[10] D. J. Watts, Small Worlds: The Dynamics of Networks 
between Order and Applications, (Princeton University 
Press, Princeton, New Jersey, 1999). 

[11] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

[12] A.-L. Barabasi, R. Albert, and H. Jeong, Physica A 281, 
69 (2000). 

[13] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002). 

[14] R. Albert, H. Jeong, and A.-L. Barabasi, Nature 406, 
378 (2000). 

[15] P. Crucitti, V. Latora, M. Marchiori, and A. Rapisarda, 

Physica A 320, 622 (2003). 
[16] P. Bak and K. Sneppen, Phys. Rev. Lett. 71, 4083 (1993). 
[17] M. Paczuski, S. Maslov, and P. Bak, Phys. Rev. E 53, 

414 (1996). 



